BELLE Preprint 2007-51 
KEK Preprint 2007-69 



o 



Measurement of Time-Dependent CP- Violating Parameters in B° — > K^Kg decays 

Y. Nakahama, 41 K. Sumisawa, 7 I. Adachi, 7 H. Aihara, 41 T. Aushev, 17 - 12 A. M. Bakich, 38 V. Balagura, 12 
E. Barbcrio, 20 I. Bcdny, 1 K. Bclous, 11 U. Bitenc, 13 A. Bondar, 1 A. Bozck, 26 M. Bracko, 19 ' 13 T. E. Browder, 6 
P. Chang, 25 Y. Chao, 25 A. Chen, 23 K.-F. Chen, 25 W. T. Chen, 23 B. G. Chcon, 5 R. Chistov, 12 I.-S. Cho, 46 
Y. Choi, 37 J. Dalseno, 20 M. Dash, 45 A. Drutskoy, 3 S. Eidclman, 1 N. Gabyshcv, 1 B. Golob, 18 ' 13 H. Ha, 15 J. Haba, 7 
K. Hara, 21 T. Hara, 31 K. Hayasaka, 21 M. Hazumi, 7 D. Heffernan, 31 Y. Hoshi, 40 Y. B. Hsiung, 25 H. J. Hyun, 16 
T. Iijima, 21 K. Inami, 21 A. Ishikawa, 34 H. Ishino, 42 R. Itch, 7 M. Iwasaki 41 Y. Iwasaki, 7 D. H. Kah, 16 H. Kaji, 21 
J. H. Kang, 46 N. Katayama, 7 H. Kawai, 2 T. Kawasaki, 28 H. Kichimi, 7 H. J. Kim, 16 H. O. Kim, 16 Y. J. Kim, 4 
£5 ! K - Kinoshita, 3 S. Korpar, 19 ' 13 P. Krizan, 18 - 13 P. Krokovny, 7 R. Kumar, 32 C. C. Kuo, 23 A. Kuzmin, 1 Y.-J. Kwon, 46 
J. Lcc, 36 J. S. Lee, 37 M. J. Lee, 36 S. E. Lee, 36 T. Lcsiak, 26 A. Limosani, 20 S.-W. Lin, 25 D. Livcntscv, 12 F. Mandl, 10 
Q ! S. McOnie, 38 T. Medvedeva, 12 W. Mitaroff, 10 K. Miyabayashi, 22 H. Miyata, 28 Y. Miyazaki, 21 R. Mizuk, 12 

G. R. Moloney, 20 E. Nakano, 30 M. Nakao, 7 H. Nakazawa, 23 Z. Natkaniec, 26 S. Nishida, 7 O. Nitoh, 44 T. Nozaki, 7 
S. Ogawa, 39 T. Ohshima, 21 S. Okuno, 14 S. L. Olscn, 6 ' 9 H. Ozaki, 7 P. Pakhlov, 12 G. Pakhlova, 12 C. W. Park, 37 

H. Park, 16 K. S. Park, 37 R. Pestotnik, 13 L. E. Piiloncn, 45 H. Sahoo, 6 Y. Sakai, 7 O. Schneider, 17 C. Schwanda, 10 
A. J. Schwartz, 3 R. Seidl, 8 ' 33 M. E. Sevior, 20 M. Shapkin, 11 C. P. Shcn, 9 H. Shibuya, 39 J.-G. Shiu, 25 B. Shwartz, 1 

, ,! J. B. Singh, 32 A. Somov, 3 S. Stanic, 29 M. Staric, 13 T. Sumiyoshi, 43 O. Tajima, 7 F. Takasaki, 7 N. Tamura, 28 

X ; M. Tanaka, 7 G. N. Taylor, 20 Y. Teramoto, 30 I. Tikhomirov, 12 K. Trabclsi, 7 T. Tsuboyama, 7 S. Uehara, 7 

V ' K. Ucno, 25 T. Uglov, 12 Y. Unno, 5 S. Uno, 7 Y. Ushiroda, 7 G. Varncr, 6 K. Vcrvink, 17 C. H. Wang, 24 P. Wang, 9 

X. L. Wang, 9 Y. Watanabc, 14 E. Won, 15 Y. Yamashita, 27 Y. Yusa, 45 Z. P. Zhang, 35 A. Zupanc, 13 and O. Zyukova 1 
r-| I (The Belle Collaboration) 

1 Budker Institute of Nuclear Physics, Novosibirsk 
I , 2 Chiba University, Chiba 

■ 3 University of Cincinnati, Cincinnati, Ohio 45221 

■^j- ' * The Graduate University for Advanced Studies, Hayama 

£C) , Hanyang University, Seoul 

£NJ ■ 6 University of Hawaii, Honolulu, Hawaii 96822 

1 High Energy Accelerator Research Organization (KEK), Tsukuba 
, 8 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 

I ■ 9 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

' 10 Institute of High Energy Physics, Vienna 

C^~) 11 Institute of High Energy Physics, Protvino 

'' • 12 Institute for Theoretical and Experimental Physics, Moscow 

. ' 13 J. Stefan Institute, Ljubljana 

, 1Jt Kanagawa University, Yokohama 

» I ■ 15 Korea University, Seoul 

d ' 6 Kyungpook National University, Taegu 

11 E cole Polytechnique Federale de Lausanne (EPFL), Lausanne 
18 Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 
19 University of Maribor, Maribor 
20 University of Melbourne, School of Physics, Victoria 3010 
Nagoya University, Nagoya 
22 Nara Women's University, Nara 
23 National Central University, Chung-li 
24 National United University, Miao Li 
25 Department of Physics, National Taiwan University, Taipei 
26 H. Niewodniczanski Institute of Nuclear Physics, Krakow 
21 Nippon Dental University, Niigata 
Niigata University, Niigata 
29 University of Nova Gorica, Nova Gorica 
30 Osaka City University, Osaka 
31 Osaka University, Osaka 
Panjab University, Chandigarh 
33 RIKEN BNL Research Center, Upton, New York 11973 



Typeset by REVTeX 



2 



3Jl Saga University, Saga 
35 University of Science and Technology of China, Hefei 
S6 Seoul National University, Seoul 
37 Sungkyunkwan University, Suwon 
38 University of Sydney, Sydney, New South Wales 
39 Toho University, Funabashi 
JfU Tohoku Gakuin University, Tagajo 
^'Department of Physics, University of Tokyo, Tokyo 
42 Tokyo Institute of Technology, Tokyo 
Jf3 Tokyo Metropolitan University, Tokyo 
" Tokyo University of Agriculture and Technology, Tokyo 
J,s Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

* Yonsei University, Seoul 

We report a measurement of the CP- violating parameters in B° — > K%K% decays based on a data 
sample of 657xlO e BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider. In this study, one neutral B meson is fully reconstructed in the 
B° — > K%Kq decay mode, and the flavor of the accompanying B meson is identified by its decay 
products. The CP-violating parameters are measured from the asymmetry in the distributions of 
the proper-time interval between the two B decays: S K o K a = — 0.38to'.7? (stat) ± 0.09(syst) and 
A k%k% = -0.38 ± 0.38(stat) ± 0.05(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



B meson decays proceeding via flavor-changing b — > 
dqq or sqq transitions are sensitive to new physics (NP) 
contributions affecting the internal quark loop diagrams. 
Such NP contributions can add new weak phases and 
subsequently cause deviations from the standard model 
(SM) expectations for CP- violating parameters [l[ . 

The decay B° — > KgKg is dominated by b — > dss 
transitions. The SM predicts that the CP-violating pa- 
rameters S K o K o and AjiO K o are zero in the limit that 
the top quark dominates and the contribution from in- 
ternal up and charm quark exchanges in the loop dia- 
gram is small Q. Measurements of S K o K a and A^k 

S3 S S 

are sensitive probes of NP [3j. The parameter S K a K o 
arises from interference between a mixing-induced am- 
plitude and a non-mixed decay amplitude while the pa- 
rameter -4^-0^-0 arises from CP violation in the decay 
amplitude itself. Experimentally, there are no prompt 
charged tracks from the B vertex for a B° — » KgKg de- 
cay, and hence the B decay vertex must be reconstructed 
by using Kg's that decay inside the vertex detector and a 
constraint on the beam interaction point (IP). Both CP- 
violating parameters have previously been measured by 
the BaBar collaboration; they obtained a large S K o K o 
value, albeit with a large statistical error In this 
Letter, we report a measurement of S K o K o and Ak 3 k% 
in B° — > KgKg decays using almost twice the statistics 
of the previous measurement. Our analysis is based on 
657 xlO 6 BB pairs collected with the Belle detector [H 
running at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider @. 

In the decay chain T(4S) — > BB — > (K%K%)f ta& , 
where one of the B mesons decays at time t K o K o to 
KgKg and the other decays at time ttag to a flavor spe- 



cific state /t a g that distinguishes between B° and B , the 
decay rate has a time dependence Q given by 

e -|At|/T B o 

V K o K o (M) = — [1 + q ■ {S K o K o sm(Am d At) 

s s 4:Tgo s s 

+ A k%k% cos(Am d At)}],(l) 

where At = t K o K o — £ tag , t b o is the B° lifetime, Am^ is 
the mass difference between the two B mass eigenstates, 
and q = +1 (-1) for / tag = B° (B°). 

At the KEKB, the T(45) resonance is produced with 
a Lorentz boost of /?7 = 0.425 nearly along the +z 
axis, which is defined as the direction antiparallel to 
the e + beamline. Since the B° and B° mesons are ap- 
proximately at rest in the T(45) center-of-mass system 
(cms), At can be determined from the displacement in 
z between the KgKg and /t ag decay vertices: At ~ 
( z k%k% -2 tag )/(/37c) = Az/(/3 7 c). 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber, an array of aerogel 
threshold Cherenkov counters, a barrel-like arrangement 
of time-of- flight scintillation counters, an electromagnetic 
calorimeter, which arc located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux return located outside of the coil is instrumented to 
detect K\ mesons and to identify muons. Two inner de- 
tector configurations were used. A 2.0 cm-radius beam 
pipe and a 3-layer SVD (SVD1) was used for the first 
data sample of 152xl0 6 BB pairs, while a 1.5 cm-radius 
beam pipe, a 4-layer SVD (SVD2) Q, and a small-cell 
inner drift chamber were used to record the remaining 
505 xlO 6 BB pairs. 

We reconstruct a Kg — > ir + ir~ candidate from a pair 
of oppositely charged tracks having \AM K o \ < 0.015 
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GeV/c 2 corresponding to three standard deviations (a), 
where AMjfO is the difference between their invariant 
mass and the nominal K s mass 0] . Both charged tracks 
are required to be displaced from the IP in the trans- 
verse (r-<p) plane by more than 100 fim. The angle in the 
transverse plane between the K s momentum vector and 
the direction defined by the K s vertex and the IP should 
be less than 50 mrad. In order to suppress incorrect com- 
binations of the two charged tracks, the mismatch in the 
z direction at the Kg vertex point for the two charged 
tracks is required to be less than 15 cm. 

To identify B° — > K S K S decay candidates, we use 
two kinematic variables: the energy difference AE = 
and the beam-energy constrained mass 
here E^ m * s ^ ne beam 
p™ s are the cms energy 



% ^bcam 



3™) 2 , W 

and 



M bc = v/PC 1 ^) 2 - (PI 
energy in the cms and E\ 
and momentum, respectively, of reconstructed B candi- 
dates. We select candidates satisfying \AE\ < 0.20 GeV 
and 5.20 GeV/c 2 < M hc < 5.30 GeV/c 2 . For the At 
fit described below, we use candidates in a signal region 
defined as \AE\ < 0.10 GeV and 5.27 GeV/c 2 < M hc < 
5.30 GeV/c 2 . We find that 0.2% of the selected events 
have multiple B° — > K S K S candidates. In such events, 
we choose the B° — > K S K S candidate having the small- 
est Y,(AM K o) 2 value. 

To suppress continuum e + e~ — > qq (q = u,d,s,c) 
events, we form a likelihood £ s i g (£ b kg) for signal (con- 
tinuum) events by combining a Fisher discriminant based 
on modified Fox- Wolfram moments [!(} with the proba- 
bility density function (PDF) for the cosine of the cms 
B Q flight direction with respect to the +z axis. The for- 
mer makes use of the difference in event shapes: signal 
events have a spherical topology while background events 
tend to be jet-like. We impose a requirement on the like- 



background is included in the systematic errors. 

The signal yield is extracted using a three-dimensional 
extended unbinncd maximum likelihood (UML) fit to 
AE-Mbc-Tl distributions for the selected candidate 
events. For the signal component, we model the AE 
(Mbc) shape using a sum of two Gaussians (a single 
Gaussian). A binned histogram is employed for the TZ 
distribution. The parameters of the Gaussians and TZ 
distribution are obtained using MC simulation. For the 
background component, the AE (M bc ) shape is modeled 
as a first-order polynomial (an ARGUS [l3j) function. 
The parameters of these functions are floated in the fit. 
The TZ background distribution is obtained from a data 
sample in the sideband region (M bc < 5.26 GeV/c 2 ). 
Possible correlations among AE, Af bc and 1Z are found 
to be negligible for signal events from the signal MC, and 
to be very small for continuum events from the data side- 
band. We include the effect of the small correlations in 
the latter in the systematic errors. The fit yields 58±11 
signal events among 476 B° — > K S K S candidate events 
in the signal region, where the error is statistical only. 
Figure Q] shows the projections of the AE, M bc and TZ 
distributions for the candidate events. 

We apply the B decay vertex reconstruction algorithm 



of Ref. 14 1 . The vertex position for a B° — > K S K S decay 



is obtained using the Kg — * ir + ir~ momentum vector 

and a constraint on the IP; the IP profile (a x ~ 100/im, 

(T y ~ 5/im) is smeared by the finite B° flight length in 

the plane perpendicular to the z axis. To reconstruct the 

B° — > KgKg decay position, both charged pions from 

at least one of the Kg's are required to have a sufficient 

number of hits in the SVD: at least one layer with hits on 

both the z and r-<j> sides and at least one additional layer 

with a hit on the z side among the other layers for SVD1, 

and at least two layers with hits on both sides for SVD2. 
lihood ratio TZ = £ sig /(£ sig + £ bkg ) that retains 89% of Both ^o> s are uscd for the vertex reconstruction if the 



the signal and rejects 71% of the continuum. 

The 6-flavor of the accompanying B meson is identified 
from inclusive properties of particles that are not asso- 
ciated with the reconstructed B° — > KgKg candidate. 
The tagging information is represented by two parame- 
ters: q, as defined in Eq. |T]), and r, which is an event- 
by-event Monte-Carlo-determined flavor-tagging dilution 
factor that ranges from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment (ll| . Candi- 
date events are selected to have r > 0.1, and are further 
divided into six r intervals. The wrong tag fraction w 
for each r interval and the differences Aw between B° 
and B° decays are determined using semileptonic and 
hadronic b — ■> c decay data 11 1 . 

The dominant background is continuum. We find the 
BB decay background contribution to be negligibly small 
using a large sample of GEANT-based Monte Carlo (MC) 
simulated events 12j. Thus, we take into account only 
signal and continuum events in the nominal fit. The un- 
certainty due to a possible contribution from BB decay 



four pions have a sufficient number of hits in the SVD. 
The typical vertex reconstruction efficiency with SVD1 
(SVD2) is determined to be 44% (61%) from the signal 
MC. The vertex position resolutions in the z direction 
with SVD1 (SVD2) arc 73/im (105/im) for the case where 
both Kg's are used for the vertex reconstruction, and 
141/im (172/xm) for the case where a single K s is used. 
The latter resolution is comparable to the /t ag vertex 
position resolution. The B decay vertex in the tag side is 
determined from well-reconstructed tracks that are not 
assigned to the B° — > K S K S decay. The typical vertex 
reconstruction efficiency for /tag decays is determined to 
be 93%. 

We determine S K a K o s and -4^o^o by performing an 
UML fit to the At distribution. For signal events, we 
use the At distribution of Eq. ([1]), modified to include 
the effect of incorrect flavor assignment. The distribu- 
tion Vko ^-o (At) is then convolved with the resolution 
function i? s ; g (Ai), which depends on the event-by-event 
vertex position errors (lBj : the dependence is calibrated 
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FIG. 1: (a) AE, (b) M bc and (c) TZ projections for the 
B° -> KgKs candidate events (a) with TZ > 0.6 and 
5.27 GeV/c 2 < M hc , (b) with TZ > 0.6 and \AE\ < 0.1 GeV 
and (c) in the signal region. The solid histogram and curves 
show the fit projections and the hatched areas show the back- 
ground component. The points with error bars are the data. 



using a B° — > J/ipKg data control sample, where the 
vertex positions are reconstructed using only a Kg and 
the IP profile [3] • We determine the following likelihood 
value for each event i: 

Pi = (l-/ol) J\SK%K%PK%K%^)^{.^ti-^) 

+(1 - f K %K%)V q - q WR qlI {AU - At')]d(At') 
+f {Poi(AU), (2) 



where the PDF V \{A£) is a broad Gaussian that repre- 
sents an outlier component with a small fraction f Q \ [15j . 
The fraction Jk°k° is the event-by-event signal fraction 
depending on AE, M^ c and TZ. We also take into account 
the r dependence of the signal fraction /x°_fst° ; we deter- 
mine the dependence using the signal MC and the side- 
band events for the signal and background components, 
respectively. For background events, the At distribution 
Vqq(At) is convolved with a function R qq (At), where the 
distribution V qq (At) is modeled as a sum of an expo- 
nential function and a delta function, and the function 
R qq (At) is the sum of two Gaussians. All parameters 
in V qq (At) and R q q(At) are determined from sideband 
events. We fix t b o and Am c i to their world-average val- 
ues @. To improve the statistical sensitivity to A K o K o, 
we also use candidate events having no At information, 
where both Kg's decay outside the SVD and we do not 
reconstruct B vertices; for these events, we use the PDF 
of Eq. ([2|) integrated over At. The only free parameters 
in the fit are S K o K o and A^o K o , which are determined 
by maximizing the likelihood function L = Y[ Pi , where 
the product is over all events. The fit to 476 B° -> K° S K° S 
candidate events, in which 216 candidate events have no 
At information, yields 

s k%k% = -0.38±£^(stat) ± 0.09(syst), and (3) 
Ak%k% = -0.38±0.38(stat)±0.05(syst), (4) 

where the systematic errors are described below. Fig- 
ure [5] shows the At distribution and raw asymmetry *4cp 
in each At interval, where .4 C p = (N+ - N-)/(N + + 
N-), and N_^r_\ is the number of candidate events with 
q= +1 (-1). 

The systematic error is primarily due to uncertain- 
ties in the parameters of i? S i g (At) (±0.06 on S K o K o and 
<0.01 on Ax°k° s )i an d uncertainties in the signal fraction 
fx k° (±0.04 on S K o K o and ±0.03 on A K o K o). We es- 

S3 S S S S 

timatc a systematic error (±0.04 on S K o K o and ±0.02 on 
Ak°k°) f° r uncertainties in the parameters of V qq (At) 
and Rqq(At), and the possible contribution of, and asym- 
metry in, the BB decay background. The other contri- 
butions to the systematic errors come from uncertain- 
ties in the wrong tag fraction (±0.02 on S K o K o, ±0.01 
on Axa K a), fit biases (±0.02, ±0.01), physics parame- 
ters (t b o and Arrid) (±0.01, ±0.01), the vertex recon- 
struction (±0.01, ±0.02), and the tag-side interference 
effect (<0. 01, ±0.03). Adding all these contributions 
in quadrature, we obtain systematic errors of ±0.09 for 
S K o K o and ±0.05 for ^l^o^o. 

Various validity checks for the measurement are per- 
formed. We measure a branching fraction for B° — > 
K°K° of [1.1 ± 0.2(stat)]xl0-^which is consistent 
with our previous measurement [171 ]. The B° lifetime 
for the B° — » KgKg candidate events is measured 
to be 1.58±0.44 ps, consistent with the world average 
value .9J. We also fit to the sideband events of the 
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the other measurement |4|. 
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FIG. 2: (a) At distribution and (b) raw asymmetry Acp for 
the B° — » Kg Kg candidate events in the signal region with 
good tags (r > 0.5) and 1Z > 0.6. In (a), the solid and dashed 
curves show the fit results with q = ±1, respectively. The 
dotted curve shows the background component. In (b), the 
solid curve shows the fit projection. 



B° — > KgKg data sample and find no CP asymme- 
try. Using MC pseudo-experiments, we find that the 
statistical errors obtained in our measurement are con- 
sistent with expectations. We apply the same procedure 
to the B° — > J / ipKg data sample without using the J/tp 
daughter tracks for the vertex reconstruction. We ob- 
tain Sj^, K o = 0.68 ± 0.06(stat), which is in agreement 
with the world average for sin2<^i 0]. We conclude that 
the vertex resolution for B° — > K^Kg decays is well- 
understood. We reconstruct 1993±53 B+ -> K° s ir + [HI 
events and, without using the charged pion of the B de- 
cay for vertex reconstruction, apply the same fit pro- 
cedure. We obtain S K a^+ = —0.13 ± 0.13(stat) and 
A K a^+ = 0.01 ± 0.06(stat), which are consistent with 
no CP asymmetry. 

In summary, we measure time-dependent CP-violating 
parameters in B° — > K^Kg decays, which are dominated 
by flavor-changing b — ► dss penguin transitions, based on 
a data sample of 657xl0 6 BB pairs recorded with the 
Belle detector. We obtain 5^^ = -0.38±^(stat) ± 
0.09(syst) and A K o K o = -0.38 ± 0.38(stat) ± 0.05(syst). 
No CP asymmetry is found for these decays. These re- 
sults are consistent with the SM prediction and also with 
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